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Abstract 

The propagation of gravitational waves or tensor perturbations in a 
perturbed Friedmann - Robertson - Walker universe filled with a perfect 
fluid is re - examined, ft is shown that while the shear and magnetic part of 
the Weyl tensor satisfy linear, homogeneous second order wave equations, 
for perfect fluids with a 7 -law equation of state satisfying | < 7 < 2, the 
electric part of the Weyl tensor satisfies a linear homogeneous third order 
equation. Solutions to these equations are obtained for a flat Friedmann - 
Robertson - Walker background and we discuss implications of this result. 

1 Introduction 

Both the covariant definition and behaviour of gravitational waves and the dual 
question of the physical interpretation of the electric (E a b) and magnetic (H a b) 
parts of the Weyl tensor are subjects of debate at present. In particular the 
question arises whether it is possible to neglect the magnetic part of the Weyl 
tensor in the Newtonian limit, and whether this is equivalent to neglecting the 
gravitational wave contribution to cosmological perturbations. It is the aim 
in this paper to extend this debate by comparing the analysis of gravitational 
waves, set in a cosmological context, using the covariant approach on the one 
hand and the more standard methods based on metric perturbations on the 
other. 

In the standard approach a second order propagation equation for the am- 
plitude of the tensor perturbations is derived either via the variation of the 
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Einstein - Hilbert action after expansion to second order, or by directly lineariz- 
ing the Einstein field equations. For example in the Bardeen formalism jjj 
the following equation of motion is obtained for the first order gauge - invariant 

(2) 

amplitude of the tensor perturbation H T , in the absence of anisotropic stress 
perturbations: [] 

+ |4 2 > + (fc 2 + 2K)H ( T 2) = , (1) 

where £ is the cosmological scale factor, K is the spatial curvature of the back- 
ground Friedmann - Robertson - Walker (FRW) model, and k is the physical 
wavenumber when K — ]T^ | . The full metric tensor perturbation is Q$ , 
where ^„ are eigenfunctions (polarization tensors) of the tensor Helmholtz 
equation on the background spatial sections. These eigenfunctions have only two 
degrees of freedom after imposing the transverse (Q^'^-v — 0) and traceless 
(Q(2)^ m = 0) conditions. 

By contrast, the covariant approach focuses on either the electric or mag- 
netic parts of the Weyl tensor O, 0], and their propagation equations when 
linearized about a FRW model (which include waves propagating in a flat - space 
background as a limiting case). In this paper we will present a rather unusual 
result: while in a vacuum both E a b and H a b satisfy second order wave equations, 
when the spacetime is non - empty, and the matter is described by a barotropic 
perfect fluid with the equation of state p = (7 — l)fx satisfying | < 7 < 2, we 
obtain a third order equation for E a b but a second order equations for both H a b 
and the shear tensor <j a t>. This is somewhat like considering Maxwell's equations 
with sources [^J, except in the gravitational case, the fundamental difference is 
that the source can never be turned off if there is matter present. 

This paper is organized as follows: in section 2 we outline the covariant char- 
acterization of gravitational waves. In section 3 we present differential equations 
in closed form which describe the propagation of gravitational waves on a FRW 
background with a perfect fluid source. Solutions to these equations for a flat 
background are given in section 4 and the paper ends with a discussion of the 
cosmological context of our results. 



1.1 Notation and Conventions 

Conventions are taken to be the same as in ^ and we take c = 1. It is assumed 
that Einstein's equations take the form G a b — ^T a b where G a b is the Einstein 
tensor, k — 8nG is the gravitational constant and T a b is the energy - momentum 
tensor of the matter, which is taken to be a perfect fluid. We use the stan- 
dard kinematical decomposition of the first covariant derivative of u a Q |||, a 
"dot" denotes the covariant derivative along the fluid flow lines and ( 3 ) V a corre- 
sponds to the 3 - dimensional covariant derivative defined by totally projecting 

lr The equation for tensor metric perturbations was first derived by Lifshitz in 1946 1141 . 
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the covariant derivative orthogonal to u a , so for example 

Tab = T a b- C u c , ^V a Xkc = h s a h f h u c T tu - s , (2) 

where 

hab = gab + U a U b (3) 

is the usual 3 - dimensional projection tensor orthogonal to u a . 

A Robertson- Walker background geometry is characterized by the shear, 
vorticity, and acceleration, together with their covariant derivatives, being at 
most first order, while the energy density and expansion are zero order, and the 
pressure can be so also |J . It follows that the electric and magnetic parts of the 
Weyl tensor: 

Eab = C ac bdU C U d H ab = \C 'acstV^ ' bdU C ' U d (4) 

are also at most first order. We will assume that in a cosmological context, 

(n + p)e > o. 

2 The covariant approach to gravitational waves 

In the covariant approach, the study of tensor perturbations was first considered 
by Hawking [jn]. In his paper, the electric part of the Weyl tensor was used as 
the variable to characterize them. Later it was suggested that the Magnetic part 
of the Weyl tensor is a better choice because it has no analogue in Newtonian 
theory where gravity is propagated instantaneously. Hence the Magnetic part 
obviously plays an important role in describing gravitational waves, but given 
the correspondence with electromagnetism [§), where neither the electric nor 
magnetic fields provide a complete description of EM waves, we suggest that 
both electric and magnetic parts of the Weyl tensor are required for a full 
understanding of tensor perturbations. Indeed it is their curls that characterize 
gravitational waves, as we will see below. 

The fully nonlinear evolution equations for E a t and H a b, with the matter 
source described by a perfect fluid, are given by (|] : 

h m a h t c E ac + h {m ail r)tsd u r H a s . 4 - 2H {t q 7 1 m)hpq UbU p + OE mt 

+h mt (a ab E ab ) - 3£ ( "V* )s - E (m s tu t)s = + p)a tm (5) 

and 

h ma h tc H ac - hJ m V r ^ d u r E a s;d + 2E^ V m)bpq u b u p 

+h mt (cr ab H ab ) + QH mt - ZH {m s o- t)s - H^ rn s uj t)s = . (6) 

Notice that just as in the electromagnetic case, the only difference between these 
equations is in the sign of the second and third terms, and the source term (here, 
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the shear) coupled to the energy density and pressure in the E a b equation. Once 
linearized about a FRW background, these equations become: 

E ab + @E ab + {3) V e H d {a rib)cdeU c + + p)a ab = (7) 

and 

Hab + OH ab - i3) V e E d {aVb)cde u c = . (8) 

On taking the time-derivative of the first and substituting from the second, it 
is the curl E and curl H terms that give rise to traveling gravitational waves, 
in analogy with the propagation of electromagnetic waves. 

It is immediately clear from the above that (given the FRW background 
evolution, which determines the zero-order coefficients in the equations), these 
equations by themselves do not close up in general, because of the last term in 
. One has to add to them the shear evolution equation, which in linearized 
form is 

Oab = (3) V( Q M b) - \h ab U c . c - |9(T ah - E ab ■ (9) 

This is crucial in leading to the third -order equation for E a b, on taking the 
time derivative of (0). 

One must also be aware of the constraint equations that have to be satisfied. 
In particular, after linearization we find the 'div E' and 'div H' equations 

E ab , b = \nX a (10) 

H ab , b = (M + PK, (11) 

where X a = /i ab /i,& is the spatial gradient of the energy density || and to a is 
the vorticity vector, which are the remaining 'Maxwell-like' equations for E a b 
and H ab - 



3 Closed evolution equations for linear gravita- 
tional waves 

In linear perturbations of FRW models, pure tensor perturbations are char- 
acterized by requiring that the following two quantities vanish to first order 




X a = => E ab - ib = 0, (12) 
u a = => H ab ;b = , (13) 

the first one excluding scalar (density) perturbations and the second one, vector 
(rotational) perturbations. The conditions that the terms on the right hand 
side vanish, following from (|l(|[ll]), are analogous to the transverse condition on 
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tensor perturbations in the metric approach, In addition, we notice that since 
the Weyl tensor is the trace - free part of the Riemann tensor, both E a b and H a b 
are trace -free, again like the tensor perturbations of the Bardeen approach. 

Given the assumed equation of state, these conditions also imply that the 
spatial gradients of the pressure and expansion and the acceleration vanish: 



Y„ 



h b 



aP,b = , Z a = h a Q, 







o, 



(14) 



the first following from the equation of state, the second from the evolution 
equations for X a (see ||), and the third from the momentum conservation 
equations. Together with ( |10|Jll|) these characterize tensor perturbations. 

Once conditions ( [l2|) and (|13|) have been imposed, one might suspect that 
both E a b and H a b would show the usual symmetry, i.e. that their propagation 
equations are the same under the substitution E a b <^> H a b- However as we will 
see below, this is not the case in general. In fact, we will demonstrate that it is 
only true when the equation of state of the background spacetime satisfies very 
special conditions. In general, as mentioned in the introduction, the equations 
for E a b and H a b are not even of the same order, the former satisfies a third 
order equation while the latter satisfies a second order one. In order to see 
this, let us begin by considering the linearized second - order equation for tensor 
perturbations, obtained by taking the time derivative of ([?]) and substituting 
from (§): 



AE ab + l@E ab + (|6 2 - 2p) E ab - |/n (37 



1 7 



Oab = 



(15) 



where AE a b — E a b — ^V 2 -E a 6- By differentiating this equation, using the 
linearised shear evolution equation (0) specialized to tensor perturbations: 



dab 



7©cr a f> — E a 



(16) 



to eliminate & a b, and using (|15| ) again to eliminate the term in o~ a b, we can elimi- 
nate the remaining shear dependence, however this still doesn't give an equation 
that is closed. This situation can be rectified if we perform a harmonic decom- 
position, expanding E a b in terms of eigenfunctions of the Laplace- Beltrami 
operator of the background FRW model: E a b = E(k)Q ab- In this way we 
can replace the Laplacian term ( 3 )\7 2 in ( |T5| ) by k 2 /£ 2 where k is the physical 
wave number if the background is flat [Q. This gives us the closed, third order 
evolution equation for the harmonic components of E a b'- 



E 



(k) 



36 



E, 



(fe) 



ie 2 



I {n + 3p) + A 



7 SB 
3 B 



E, 



(k) 



where 



A 



AB 



w ^lQA-B 



E, 



(*0 



0, 



a= |e 2 -2 P - 



(17) 
(18) 
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and 

B = -\m [(3 7 -2)0- . (19) 

In contrast, the propagation equation for H ab is relatively simple: 

AH ab + lQH ab + (|9 2 - 2p) H ab = 0. (20) 

The reason why the equation for H ab is second order is due to the constraint 
equation || (with vanishing vorticity): 

H ab = - (3) VV d (a%)cde u c . (21) 

This allows the shear term - arising on taking the time derivative of (^) and 
substituting from (R) - to be replaced in terms of H ab . In the case of the 
equation for E ab (|15|), the shear term cannot be removed without differentiating 
again, leading to the higher order equation. 

The shear also satisfies a second order equation. This follows by differenti- 
ating equation ( |l6| ) and using (0) and (^lj) to substitute for E ab and H ab : 

Aa ab + I Qa ab + [±9 2 + \p (9 7 - 1)] a ab = . (22) 

By comparing the equations ([l5|) for E a b and (|20|) for H ab , it is clear that 
they will be identical in form iff the equation of state (assumed to satisfy 7^ 7^ 0) 
satisfies 

7 -i7(37-2) = 0, (23) 

for then the last term in ( |l5| ) will vanish. In the case of a barotropic perfect 
fluid where 7 is a constant, two possibilities exist: either j — 0^-p+p = 
corresponding to a false vacuum, already rejected, or 7 = | which represents 
the coasting solution and is also non-physical. If 7 is allowed to vary, then ( f23|) 
will also be satisfied provided 

7 = 2 2 . , (24) 

P — eXp 

where /3 is an arbitrary constant. However this will usually have to be rejected 
because it will then, for early or late times, exceed the standard bounds < 
dp/d[i < 1 demanded for physical reasons H (note that when 7 varies, it is 
no longer true that dp/dfi — 7). For normal matter satisfying | < 7 < 2, the 
equations for E ab and H ab will be very different. 

The appearance of a third order equation for the electric part of the Weyl 
tensor E ab is very surprising. First of all the standard theory discussed in sec- 
tion ^ gives a second order evolution equation (Q) , and secondly, force laws are 
expected to be formulated as second order evolution equations. We consider this 
again in section 6. Hawking's covariant analysis |llj gave a second order equa- 
tion, because he assumed a vacuum condition, instead of consistently sticking 
to a cosmological context (as assumed here). We consider this again in section 
5.1. 
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4 Solutions for a flat FRW background 



At first glance there seems to be an inconsistency with waves mediated by E a t, 
and -ff a b, where one is governed by a second order equation and the other by a 
third order equation. However consistency of the equations has been carefully 
checked |], We need then to explore the nature of the solutions, and also 
check whether there is some hidden symmetry in the equations that makes them 
compatible. 

For simplicity we will solve the above equations only in the case of a flat 



FRW background (K = 0). Using the conformal time variable rj: 



dri 1 

dt 



the 



equations for the harmonic components of the shear and magnetic part of the 
Weyl tensor become 



f^)~[^ 2 (3M-4)-fe 2 ] a w =0, 



(25) 



and 



H" k ) 



■ H [k) 



for a flat background the scale factor 
given by: 



[2^ 2 ( 7 - 2) - fc 2 ] H W = . (26) 
energy density /i and expansion O are 



67 



I OC ?7 3 T-2 5 fi (XT) 3 T~ 2 , 

Substituting the background parameters into 
gives the following general solutions 



37 



CX 77 3 7-2 

(H) and (H) 



(27) 

and integrating, 



37-10 r 



J (k) 



^2(37-2) I J v (kTl)C[l] 



(28) 



and 



37-14 



[Ju(krj)D^ + N v {krj)D^ 



(29) 



where C and D are arbitrary constants and J v , N v denote the Bessel functions 
of the first and second kind. 

For large scales (fc — > 0), we obtain power- law solutions: 



T{k) 



22^4 6_ 



and 



y (fc) 

3 (7-2) 



(fc) 



(30) 



H {k) = D^rj ^ +D$r, 37-2 . (31) 
In the case of a flat background, the third order equation for Eha ( [l7|) becomes: 



E(k) + (4 + 3 7 )^ ( " fe) + [(4 7 + 6) / i£ 2 + 



+ {i[16- 7 (37-2)] M £ 2 + | 7 fc 2 }£9£; i 



H*0 



(32) 
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This time, a general solution could not be found, however for large scale modes 
we obtain: 



E {k) = r, 3,-2 37 - 2 E m + r,37-a B gj . (33) 

The first two of these modes were previously known || [l(| . The third one, due 
to the third order nature of the equation for En.), has to be unphysical as it is 
not a solution of the original system of first order equations. The point is an 
essential one: one cannot satisfy one or other of these equations in isolation, 
one has to solve the set as a whole. In this case, we can solve first for <T a j, (a 
second order equation), as above, and then determine E a f, from ( |l6| ) by simply 
differentiating this solution. This will determine solutions to the third order 
equation for E a b] but there will be two modes, not three: 

3(7+2) 4 
E W = V - 37-2 E§+t,-*i-*E<». (34) 

This shows that in fact the effective equation governing the solution for E a b, 
considered as a solution to the set of equations as a whole, is second order, 
corresponding to the previously known solutions, even though one apparently 
cannot write down a separate second order equation for E a b alone. 



5 The cosmological context 

We have emphasized here that we are considering the problem of gravitational 
waves in a cosmological setting. 



5.1 Averaging 

The key point in the analysis has been a non-vacuum assumption. Is this real- 
istic? 

As discussed above, the evolution equations for E a b and H a b are the same 
when fx + p = 0. Now for astrophysical sources of gravitational waves, such as 
pulsars or merging binary systems in our own galaxy, the vacuum assumption is 
likely to be good for much of the region the wave travels in, since the intergalactic 
and interstellar media are relatively tenuous. Waves with wavelength less than 
1 Mpc will in fact travel through a vacuum for much of the time; and those of 
wavelengths less than 1 pc will be in an effective vacuum for almost all the time. 
Thus in this case the Hawking analysis will be adequate; both E a b and H a b will 
satisfy the same second-order wave equation. 

However cosmological scale gravitational waves, and certainly those with 
wavelengths of the order of the horizon size, will experience a non- vacuum 
spacetime background all the time. Thus our problem is one of averaging: on 
what scale does the geometry of the universe start to approximate that of a 
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Robertson- Walker universe? M On smaller scales, the vacuum approximation 
will be acceptable for much but not all the time; on larger scales it is not, and 
the third order equation discussed above for E a b applies. 

In particular, we can ask the following question: What happens to E a b 
for small-scale waves if, after traveling in vacuum, it comes across a non-zero 
matter concentration, for example in the case of a gravitational wave crossing a 
galaxy cluster? The wave equation for the waves will change then from second 
to third order. However perhaps this is not too drastic, as the effective order of 
the set of equations will remain second order, as discussed above. Nevertheless 
the junction conditions at such a change need careful consideration, because 
the standard assumptions allow a discontinuity of the curvature tensor at a 
boundary of this kind. We will not pursue this further here, but note it as 
worthy of investigation. 

5.2 Sharp phase transitions 

The same kind of issue arises in the early universe, in a different context. 

If the early universe went through an inflationary phase, then there must 
have been a transition at the end of inflation from the de - Sitter phase to a ra- 
diation dominated phase, often treated as an instantaneous 3- surface, although 
in reality reheating is not instantaneous. The question we pose is as follows: 
can there be amplification of perturbations across this boundary? The standard 
answer is that density perturbations can be amplified , but that tensor pertur- 
bations cannot. This is simply a result of the Darmois junction conditions which 
require that the three - metric and the extrinsic curvature of the three - metric 
be continuous across the boundary. 

The key point for us is that for tensor perturbations, this theorem of no- 
amplification does not carry through in the covariant case since both the electric 
and magnetic parts of the Weyl tensor will depend on the second derivatives of 
the metric perturbations, which are unconstrained by the junction conditions. 
Thus we cannot agree with the conclusion that gravitational waves cannot be 
amplified across a sharp phase transition. We will not discuss this issue further 
here; it is raised because it is similar to the situation just discussed. 

6 Conclusion 

It is interesting that the evolution equations for E a b and H a b are so different, 
because this implies a radical break with the analogy of source-free electromag- 
netism, which had been thought to carry through almost completely, at least 
qualitatively. This is assuming of course that the evolution equation for E can- 
not be reduced to second order (the 'hidden symmetry' option). We have not 
found a way to do so; the third order equation appears to us to be genuinely 
third order. 



Clearly it is interesting to locate the conditions when the evolution equations 
for E a b and H a b reduce to the same second order equation. We found that this 
happens (i) when (/i+p) = 0, i.e. in vacuum and in de Sitter spacetime, (ii) for 
the exceptional equation of state p = —(1/3)//, (iii) and the unphysical variation 
of 7 given by (]24|). We exclude all these cases in a realistic cosmological context, 
for long- wavelength waves. 

Our analysis shows that for large scale gravitational waves in a universe 
model with realistic matter, E a b and H a f, obey quite different equations, yet 
these are consistent with each other. This is like the situation for Maxwell's 
equations with a source, as is shown in the accompanying paper Although 
the wave equation for E a b is third order, it has the same characteristics as the 
second order equation for H a b, and both can be solved from a single 'potential', 
namely the shear (obeying a second order wave equation). This presumably is 
how the Bardeen analysis leads only to a second order equation: that equation 
is also for a potential, and the equation satisfied by the Weyl tensor itself is not 
written down. 
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